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Abstract. We present results of mid-infrared A = 5.0 — 16.5 fim spectrophotometric imaging of the starburst 
galaxies M82, NGC253, and NGC 1808 from the ISOCAM instrument on board the Infrared Space Observatory. 
The mid- infrared spectra of the three galaxies are very similar in terms of features present. The A ^ 11 fim 
continuum attributed to very small dust grains (VSGs) exhibits a large spread in intensity relative to the short- 
wavelength emission. We find that the 15 /^m dust continuum flux density correlates well with the fine-structure 
[Aril] 6.99 /im line fiux and thus provides a good quantitative indicator of the level of star formation activity. 
By contrast, the A = 5 — 11 /im region dominated by emission from polycyclic aromatic hydrocarbons (PAHs) 
has a nearly invariant shape. Variations in the relative intensities of the PAH features are nevertheless observed, 
at the 20% — 100% level. We illustrate extinction effects on the shape of the mid-infrared spectrum of obscured 
starbursts, emphasizing the differences depending on the applicable extinction law and the consequences for the 
interpretation of PAH ratios and extinction estimates. The relative spatial distributions of the PAH, VSG, and 
[Ar ll] 6.99 fim emission between the three galaxies exhibit remarkable differences. The ^ 1 kpc size of the mid- 
infrared source is much smaller than the optical extent of our sample galaxies and 70% — 100% of the IRAS 12 fim 
flux is recovered within the ISOCAM < 1.5 arcmin'^ field of view, indicating that the nuclear starburst dominates 
the total mid-infrared emission while diffuse light from quiescent disk star formation contributes little. 
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1. Introduction 

As hosts of vigorous episodes of star formation activity, 
starburst galaxies in the local Universe offer a unique 
opportunity to study the star formation process in ex- 
treme environments reminiscent of the conditions presum- 
ably prevailing in primordial galaxies. Studies in recent 
years have underscored the prominent role of starbursts 
in galaxy formation and evolution as well as their signifi- 
cant contribution to the extragalactic background, and to 
the chemical enrichment and heating of the intergalactic 



medium up to the highest redshifts (e.g. Steidel et al. 

|l996t Small et al. |1997|; Franx et al. 



1996 



1998; Elbaz et al. 2002). Understanding the 
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starburst phenomenon is thus a key issue of local and cos- 
mological relevance. 

Nuclear starburst galaxies generally emit the bulk of 
their luminosity in the infrared, primarily as reprocessed 
radiation by interstellar dust grains heated by important 
populations of hot massive stars. In addition, a substan- 
tial fraction - if not most - of the star formation in star- 
burst systems is heavily obscured at optical and ultraviolet 
wavelengths by large amounts of dust. This has been par- 
ticularly dramatically illustrated by mid-infrared imaging 



of the colliding pair NGC 4038/4039 (Mirabel et al. |T99§). 
Studies at infrared wavelengths are thus crucial in defin- 
ing the properties of starburst galaxies and the interplay 
between the primary energy sources and the reprocessing 
material. The growing evidence for the presence and sig- 
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nificant role of dust in distant galaxies further emphasizes 
the need for infrared investigations of nearby dust-rich 
templates (e.g. Rowan- Robinson 2001; Blain et al. 2002; 
and references therein). However, except for a limited sam- 
ple, the detailed spatial and spectral energy distributions 
of starburst galaxies in this regime are still poorly known. 
Progress has been hindered by instrumental limitations, 
and by the difficulties inherent to ground-based and air- 
borne observations at these wavelengths. 

The Infrared Space Observatory {ISO; Kessler et 
al. 1996] ) has revolutionized the field of infrared as- 
tronomy (see the reviews by Cesarsky & Sauvage 1999; 
Genzel & Cesarsky 200C), with unprecedented capabili- 
ties compared to its predecessor the Infrared Astronomical 
Satellite (IRAS). In the mid-infrared (MIR, A = 5.0 - 



16.5 /xm), the camera ISOCAM (Cesarsky et al. 1996) pro- 
vided an improvement in sensitivity and spatial resolution 
by factors of ^ 1000 and ~ 60, respectively. In addition 
to broad- and narrow-band filters, ISOCAM featured con- 
tinuously variable filters (CVFs) allowing arcsecond-scale 
spectrophotometric imaging with full wavelength cover- 
age at a resolution oi R = A/AA ~ 40. The MIR regime 
contains a variety of key features tracing different compo- 
nents of the interstellar medium (ISM) and is subject to 
relatively little extinction by dust (As^iy^m < 0.1 Ay). 
Spectrophotometric observations permit studies of the 
most prominent features to characterize the spatial dis- 
tribution and infer the properties of the ISM and exciting 
sources, deep into obscured star-forming regions. 

In this paper, we present results of ISOCAM CVF ob- 
servations of M82, NGC 253, and NGC 1808. M82 and 
NGC 253 are two of the nearest and brightest starburst 
galaxies (at distances of £> = 3.3 and 2.5 Mpc, respec- 
tively, 1" = 16 and 12 p c; Freedman & Madore |1988 ; 
Davidge & Pritchet 1990). Both have been extensively 
studied in all accessible spectral regions and have often 
been considered as prototypical objects for the starburst 
phenomenon (see Telesco 1988, Rieke et al. 1993, and 



Engelbracht et al. |1998 for reviews). NGC 1808 is a more 
distant system {D = 10.9 Mpc for Hp = 75 kms"^ Mpc"\ 
1" = 53 pc; Sandage & Tammann 1987) which closely re- 



sembles M 82 a nd N GC 253 in its global properties (e.g. 

1995|) . In all three 



Dahlem et al. 1990; Junkes et al 



galaxies, starburst activity takes place within the central 
0.1 — 1 kpc regions which are severely obscured at opti- 
cal and ultraviolet wavelengths. The large extinction lev- 
els likely result in part from high disk inclination angles 
of 80°, 78°, and 57° for M 82 , NGC 253, and NGC 1808 
(Gotz et al. |199G| ; Pence |l98l[ Reif et al. |l982D . 

At comparable infrared luminosities^ Lir = 
is-ioooMm of 3.7 X 101°, 1.4 X 10^°, and 3.4 x 10^° Lq for 
M82, NGC 253, and NGC 1808, these galaxies also share 
key morphological features including a kiloparsec-scale 
stellar bar, large amounts of molecular gas in the central 
~ 1 kpc predominantly concentrated in circumnuclear 



ring-, spiral-, or bar-like structures, populations of 
compact and luminous "super star clusters" as well 
as of compact non-thermal radio sources associated 
with supernova remnants, and a large-scale outflowing 
starburst wind (Forster Schreiber et al. 2001; Engelbracht 
et al. 1998; Tacconi-Garman et al. 1996; and references 



therein). The triggering of starburst activity in M82 
is generally attributed to gravitational interaction as 
evidenced by the extended filamentary tidal features in 
Hi gas distribution threading M82 and its neighbours 
M81 and NGC 3077 or, alternatively, to the stellar bar 
which may itself have been induced by the interaction 
(e.g. Yun et al. |l99|, [l994| ). In NGC 253, the bar itself 
appears to be the primary mechanism responsible for 
the onset of starburst activity as there are no obvious 
signatures of interaction with a nearby companion, 
although a past merger or accretion event involving a 
small galaxy has been suggested based on kinematical 
evi dence (e.g. Anantharamaiah & Goss 1996; Boker et 
al. |1998[ ). Tidal interaction with NGC 1792 has been 
proposed for NGC 18 08 ba sed on circumstant ial ev idence 
(e.g. Dahlem et al. |l990t Koribalski et al. |1993| ). The 
recent high sensitivity Hi observations by Dahlem et al. 
(2001) do not indicate the presence of any tidal feature, 
dwarf galaxy or other debris that could support a close 
passage in the past, leaving the stellar bar as more 
probable trigger or a possible merger/accretion event. 

In view of their similar nature, we focus on the com- 
parison of the spatial and spectral properties of M82, 
NGC 253, and NGC 1808 at MIR wavelengths. Though 
small, this sample also allows us to probe the transition 
regime between normal spiral and irregular galaxies and 
the more extreme luminous and ultraluminous infrared 
galaxies (LIRGs and ULIRGs for which Lir > 10^^ and 
> 10^^ Lq, respectively). Section || describes the observa- 
tions and data reduction procedure, and Sect. |^ presents 
the results. Section^ discusses the origin and spatial disti- 
bution of the continuum and emission features. Section |^ 
addresses the issues of extinction effects, variations in 
spectral properties, and indicators of star formation ac- 
tivity. Section ^ summarizes the paper. 

2. Observations and data reduction 

We observed M82, NGC 253, and NGC 1808 with 



^ as computed from the I RAS fluxes following the prescrip- 
tion of Sanders & Mirabel (|1996|). 



ISOC AM (C esarsky et al. |1996D on board ISO (Kessler 
et al. 1996 ) as part of the guaranteed time program 
CAMACTIV (P.I. I. F. Mirabel). We used the CVFs 
to cover the entire A = 5.0 — 16.5 /im range with a 
resolution of i? = A/AA w 35 to 45 and increments 
SX — 0.06 to 0.1 fj,m from short to long wavelengths. For 
all three sources, twelve single-frame exposures of 2.1 s 
were recorded per wavelength channel, for total integra- 
tion times of about 1^05™. We selected the 3" pixel" ^ scale 
for M82 and NGC 1808, and the 1.5"pixel-i scale for 
NGC 253; total fields of view of 96" x 96" and 48" x 48" 
are thus covered by the 32 x 32 pixels detector array. 
Table |l| presents the details of the observations. The ar- 
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Fig. 1. The ISOCAM CVF field of view, indicated by a square, is superimposed on publicly available images in the 
V band for M82 (left), and R band for NGC 253 and NGC 1808 (middle and right). 

Table 1. Log of the observations 



Source 


Date 


Revolution 


Q2000 


(52000 


P.A." 


Pixel scale 


f- ^ 

tint 








. m . 


(° : ' : ") 


(degree) 


(arcsec) 


(s) 


M82 


1996 Mar 19 


123 


09:55:51.8 


+69:40:45.8 


212 


3.0 


25.2 


NGC 253 


1997 Jan 03 


414 


00:47:32.9 


-25:17:18.3 


290 


1.5 


25.2 


NGC 1808 


1997 Oct 12 


697 


05:07:42.3 


-37:30:46.1 


64 


3.0 


25.2 



" Position angle of the vertical axis of the ISOCAM array relative to the north and increasing counterclockwise. 
Total integration time per wavelength channel. 



eas observed with ISOCAM cover almost entirely the MIR 
sources in each galaxy (Sect. ^) and are indicated on opti- 
cal images in Fig. |^. This emphasizes the small MIR source 
sizes compared to the optical extent of the galaxies and 
that most of the MIR emission originates in their nuclear 
regions. 

We reduced the data using the CAM Interactive 
Analysis package (CIA)|^ following the procedures de- 
scribed in the ISOCAM User's Manual (Delaney |1997D . 
We first subtracted the dark current using a model of the 
sec ular e volution of ISOCAM's dark current (Biviano et 
al. 1997), removed cosmic ray hits by applying a multires- 
olution median filtering technique (Starck et al. 1999| ), 
and corrected for transient detector memory effects us- 
ing t he "Fouks-Schubert" algorithm (Coulais & Abergel 
200C ) . We then combined the individual exposures for each 
wavelength channel with the help of a shift-and-add algo- 
rithm accounting for the satellite jitter motions (typical 
amplitude of k, 1") and performed absolute flux calibra- 
tion based on the sensitivity calibration files provided in 
CIA version 3.0. 

For proper comparison of fiuxes and maps obtained at 
different wavelengths, we smoothed the CVF channels to 
a common angular resolution and registered them appro- 
priately. The point spread function (PSF) varies throuh- 



^ CIA is a joint development by the ESA Astrophysics 
Division and the ISOCAM consortium led by the ISOCAM 
P.I. C. Cesarsky, Direction des Sciences de la Matiere, C.E.A. 
Saclay, France. 



gout the CVF scans, with theoretical full-width at half- 
maximum FWHM = 3.1" - 5.6" for the 3"pixel-i scale 
and 2.0" — 5.2" for the 1.5"pixel~^ scale, over the range 
5.5 — 16.5 fim. We convolved the data to a final resolution 
of 5.6" for M82 and NGC 1808, and 5.2" for NGC 253, 
corresponding to a linear resolution of 90, 297, and 62 pc, 
respectively. We computed the transfer functions (close 
to Gaussians in profile) using the library of theoretical 
PSFs for ISOCAM and a CLEAN algorithm. The small- 
scale structure is thus smeared out in the higher resolution 
channels but no important spatial feature is lost. 

The source position varies slightly with wavelength due 
to telescope motions and changes in optical path, with the 
largest offset occurring near 9.2 /im between the short and 
long wavelength segments of the CVFs. Since no unre- 
solved point source is detected in the data, we determined 
the position offsets from cross-correlation over the entire 
field of view, taking as references images integrated over 
various wavelength intervals. Effects due to the different 
spatial distributions of emission features were easily iden- 
tified and ignored in the determination of the instrumental 
shifts. We fitted straight lines to the cross-correlation off- 
sets for each CVF segment. For the 3"pixel^"'^ scale, the 
horizontal, or x-axis offsets relative to 9.2 /xm in the short 
and long wavelength segments did not exceed 0.3 pixel 
and 0.5 pixel, respectively, with a jump of 0.7 pixel at the 
transition wavelength. Shifts along the j/-axis were negligi- 
ble throughout the entire wavelength range (< 0.2 pixel). 
For the 1.5"pixel~^ scale, both x- and y-axis offsets were 
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< 0.3 pixel within each CVF segment, with jumps at 
9.2 fim of 0.2 pixel only. 

We did not attempt to correct for instrumental flat- 
field, stray light, and ghost images generated by multiple 
reflections between the detector and filters since appropri- 
ate calibration files for the CVF mode were not available. 
Based on theoretical models for point-like and extended 
sources as well as observations of stars and of the zodia- 



cal hght (Biviano et al. [19984 |1998b| Okumura |2000[ ), the 
combined flat-field, straylight, and ghost effects integrated 
over the entire detector array vary between « 35% and 
10% of the source flux from 5 to 16 /im, and exhibit quali- 
tatively similar spatial structures at different wavelengths. 
The r.m.s. noise estimated from existing flat-field images 
obtained with ISOCAM's broad-band filters is 10% for 



both pixel scales (Biviano et al. |l998b|) . 



Based on available calibration accuracy reports 



(Blommaert & Cesarsky 199g and references therein), the 
systematic errors of the data processing and flux cali- 
bration described above sum up to total uncertainties on 
the absolute photometry of « 50%. For our analysis, we 
will consider formal effective uncertainties on the relative 
fluxes estimated as follows. We measured the dispersion 
around the mean flux among the individual reduced expo- 
sures for each wavelength channel at each detector pixel, 
included an additional 10% r.m.s. noise for intrinsic pixel- 
to-pixel variations in sensitivity, and accounted for un- 
certainties in relative spectral calibration. For M82 and 
NGC 253, the latter were evaluated from the comparison 
with data from the ISO Short Wavelength Spectrometer 
(SWS; de Graauw et al. I99t) in Sect. ^, amounting to 
5% and 11%, respectively. For NGC 1808, very low flux 
levels and high noise around 10 /im., below 6 /xm, and out- 
side of the main emission region led to poor correction 
for transient memory effects because of the difficulty in 
determining the initial and stabilized fluxes and of the 
bad behaviour of the algorithm wherever fluxes become 
negative. In this case, we estimated the errors in relative 
spectral calibration by comparing results when applying 
the transient correction to the data as they were and to 
the data where all values lower than 0.01 were set to 0.01, 
adding a conservative 10%. The median and average dif- 
ferences are of 16% and a factor of 2. Experimentation 
showed that the choice of the threshold is of little conse- 
quence as long as it is sufficiently smaller than the signal 
where the transient correction is satisfactory. 

Overall, the formal effective uncertainties per pixel and 
wavelength channel within the brighter emission regions of 
interest (defined by the "starburst cores" in Sect. ^ have 
median values of 18%, 29%, and 45% for M82, NGC 253, 
and NGC 1808. Computed over the entire ISOCAM field 
of view, the medians increase to 32%, 88%, and a factor 
of 2.6 as a result of the higher noise towards the edges of 
the array where little flux is detected. 



3. Results 

3.1. Mid-infrared spectra and feature identification 

We extracted from the reduced ISOCAM data cubes spec- 
tra integrated over various regions: a "starburst core" cen- 
tered on each galaxy's nucleus^ and covering most of the 
emission source, and a smaller region at the position of 
the peak observed for the integrated MIR emission. For 
M82 and NGC 253, we further included a "disk" corre- 
sponding to an annulus outside of the starburst core and 
the ISOCAM field of view. Due to the higher noise level 
in the NGC 1808 data, only the starburst core and MIR 
peak spectra are of sufficient quality for analysis. The disk, 
core, and MIR peak regions were chosen to sample rela- 
tively quiescent to intense star-forming activity as traced 
by e.g. the 15 /im continuum emission (see Sects. 3.2 and 
|3.2D . The starburst cores correspond to the apertures gen- 
erally used in previous studies when referring to the global 
properties of the starburst in each galaxy. 

Figure || indicates the selected regions on maps of the 
14.8 — 15.2 /im emission. Figure || shows the spectra, all 
normalized to a total 6.0 — 6.6 /im flux density of unity 
to facilitate comparison of the relative strength of the 
emission features at A <i 11 /im and of the continuum in- 
tensity at longer wavelengths. The formal effective uncer- 
tainties are plotted along with the spectra. Table || gives 
the parameters for the synthetic apertures, the normaliz- 
ing fluxes, and the average and median uncertainties for 
each spectrum. Figure || also shows for M 82 and NGC 253 
the average spectrum of all individual resolution elements 
within the ISOCAM field of view, normalized as described 
above, together with the dispersion around the mean, the 
full range observed, and the typical (median) uncertain- 
ties. The resolution elements for M 82 and NGC 253 corre- 
spond to rebinned pixels of size 2x2 and 4 x 4 in original 
detector pixels, respectively. 

All spectra look very similar and exhibit the "classi- 
cal" characteristics observed towards star-forming regions 
and galaxies: conspicuous broad emission features (the 
"unidentified infrared bands" or UIBs), a featureless con- 
tinuum rising importantly at A ^ 11 /im, and an apparent 
dip near 10 /xm. The UIBs are attributed to a family of 
particles, the nature of which still is debated, stochasti- 
cally heated by single ultraviolet photons while the long- 
wavelength continuum is ascribed to very small dust grains 
between the transient heating and thermal equilibrium 
regimes depending on grain properties and radiation field 



intensity (e.g. Leger et al. 198£; AUamandola et al. 1989; 
Desert et al. |1990|; Duley & Williams |l99l|; Tielens et 



al. 1999; see also reviews by Puget & Leger 198E; Cesarsky 
& Sauvage |1999| Genzel & Cesarsky ^OOOl ). We will here- 
after refer to these ISM components as "PAHs," adopt- 
ing the currently popular model in which the UIB carriers 
consist of polycyclic aromatic hydrocarbon molecules, and 
"VSGs." 



as defined by the peak of the near-infrared 2 /im emission; 
see table for coordinates and references. 
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Fig. 2. Selected regions in M 82 (top), NGC 253 (middle), 
and NGC 1808 (bottom), shown on A = 14.8 - 15.2 /im 
continuum maps. The starburst core and MIR peak re- 
gions are enclosed within the black and white circles, re- 
spectively. The disk regions for M 82 and NGC 253 corre- 
spond to the hatched areas. The parameters for all regions 
are given in table ^. The horizontal bar and the filled circle 
at the bottom of each panel indicate the physical linear 
scale and the FWHM of the PSF. 



A number of weaker emission features are also detected 
in the spectra of Fig. ^ However, their identification is 
problematic at such low spectral resolution due to pos- 
sible feature blends. To emphasize this point and secure 
the identifications, Fig. ^ shows the spectra of M 82 and 
NGC 253 obtained at R - 500 - 1000 with the ISO-SWS 
(from Sturm et al. pOOO and Forster Schreiber et al. 2001). 
The SWS data are also compared with ISOCAM spectra 



taken in the same apertures, after convolution to the same 
spectral resolution. The ISOCAM and SWS data agree 
very well, confirming the accuracy of the absolute and 
relative flux calibration for both instruments. For M82, 
the differences are 5% on average (24% at most). For 
NGC 253, they are of 11% on average (30% at most) at 
A > 5.4 /im while they reach a factor of two at shorter 
wavelengths, probably due to residual transient memory 
effects. 



From this comparison (see also Sturm et al. 2000), 
the features at 6.2, 7.7, 8.6, and 11.3 /im are unambigu- 
ously identified with PAH emission. The broad feature 
near 12.7 fj,m clearly results from the blending of the PAH 
12.7 /im band and of the [Neii] 12.81 /tm fine-structure 
line. The [Nelll] 15.56 /tm line is also blended with the 
nearby PAH 15.7 //m feature. From the SWS data, and 
with the continuum and integration bandpasses given be- 
low, the PAH 12.7 /im accounts for about 50% of the flux 
in the 12.7 /im blend for both galaxies. The PAH 15.7 /tm 
contains about 30% of the flux in the 15.6 /im blend for 
M 82 (the noisier SWS spectrum of NGC 253 makes an 
estimate difficult). The least contaminated fine-structure 
line detected with ISOCAM is [Ar ii] 6.99 /im, with the un- 
derlying PAH at 7.0 nm and the H2 — S{5) rotational 
line at 6.91 /im contributing « 20% and 5%, respectively, 
to the blend flux in both galaxies. The weak features at 
5.65, 13.55, and 14.25 /tm are identified with PAH bands; 
the latter is definitely not due to the high excitation [Ne v] 
14.32 /im line which is undetected in the SWS data. The 
final identifications in the ISOCAM spectra are given in 
Fig. ^. Because of the excellent correspondence between 
features seen in all three galaxies, the identifications for 
M 82 and NGC 253 are adopted for NGC 1808 as well. 

We performed various continuum and feature measure- 
ments on the spectra. Table K[gives the intervals used in 
the computations and tables^ and ^ report the results. 
The flux densities in the 6.0 — 9.0 /tm and 13.5 — 15.0 /tm 
bands are dominated by emission from PAHs and VSGs, 
respectively. We chose the narrow continuum bands cen- 
tered at 5.5 and 15.0 /im so as to minimize the contribution 
from PAHs and other emission lines. In particular, we ob- 
tain the same 5.5 /im flux densities within « 5% from the 
ISOCAM data and from the higher resolution SWS spec- 
tra for the SWS field of view in both M 82 and NGC 253, 
indicating negligible contribution from the adjacent PAH 
5.65 /im feature and [Fe 11] 5.34 /im line in the lower resolu- 
tion data. We measured the fluxes in the PAH features at 
6.2, 7.7, 8.6, and 11.3 /im, in the [Am] 6.99 /im line, and 
in the PAH 12.7 /tm + [Neil] 12.81 /im and PAH 15.7 /im 
+ [Neiii] 15.56 /tm blends by integrating the flux under 
the feature profiles after subtracting a continuum baseline 
obtained by linear interpolation between adjacent spec- 
tral intervals. More sophisticated methods such as profile 



fitting (e.g. Uchida et al. 2000) are not necessary for our 
purposes and are difficult to apply to the lower signal-to- 
noise (S/N) ratio data of individual pixels in generating 
linemaps (Sect. 3.2). 
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Table 2. Synthetic apertures for selected regions 







PoQit inn 


Aa " 


AS " 


A npvtnvp 


J6. 0—6.6 


T TtiPPvtaintipQ 










(circsec) 


( avr^pp ) 

\ CAr± VjOtI'Vj J 




(Vo) 


M 82 


PAM POV 

\^iT.lVi 1: V. ./ V 




u.u 




R^ V 81 at P A of 919° 




(I) 




Disk 


Nucleus 


0.0 


0.0 


20 < r < 30 


9.37 ± 0.06 


3 (2) 






l\ 1 1 P 1 PI 1 Q 


0.0 


0.0 


r < 15 


19 5 + 02 


2 (2'] 




MIR peak 


MIR Peak 


-3.0 


-1.5 


r < 6 


5.74 + 0.10 


5 (5) 


NGC 253 


CAM FOV 




-4.0 


+1.0 


40.5 X 40.5 at P.A. 290° 


8.34 + 0.73 


20 (3) 




Disk 


Nucleus 


0.0 


0.0 


10 < r- < 15 


1.75 + 0.09 


11 (4) 




Core 


Nucleus 


0.0 


0.0 


r < 7.5 


3.90 + 0.05 


4(3) 




MIR peak 


MIR Peak 


0.0 


0.0 


r < 3 


1.14 + 0.03 


8(7) 


NGC 1808 


Core 


Nucleus 


0.0 


0.0 


r < 15 


1.77 + 0.15 


57 (4) 




MIR peak 


MIR Peak 


+1.5 


-1.5 


r < 6 


0.74 + 0.02 


9(8) 



Right ascension and declination o ffsets relative to the position of the nucleus of each galaxy. For M 82: Q2000 = 09''55™52:2, 



+69°40'46'.'6 (Dietz et al. 1986). For NGC 253: 02000 = 00'' 47'"33 n, 82 



-25°17'18'.'3 (Sams et al. 1994). For 



NGC 1808: 02000 = 05''07"'42!3, ^2000 = -37°30'46'.'3 (CoUison et al. |1994|). 
'' Integrated flux density between 6.0 and 6.6 used for the normalization of the spectra plotted in Fig. ^. The integration 
area reported for the ISOCAM field of view of M 82 and NGC 253 excludes the edges which were not illuminated or are 
highly noisy. Quoted uncertainties are computed from the formal effective uncertainties for individual pixels and wavelength 
channels (see Sect. 

Average of the formal effective uncertainties per wavelength channel for the integrated spectra (median values are given in 
parenthesis). 



We did not compute fluxes for the individual [Nell] 
12.81 /ini and [Ne ill] 15.56 /im lines. The wavelength sam- 
pling is too coarse for reliable profile decomposition, 
with the feature peaks of each blend sampled by ad- 
jacent wavelength channels. In addition, our attempts 
to subtract the PAH contribution by attributing an 
excess in the blend 12. 7/im/PAH 11.3 /im ratio to the 
[Neil] 12.81 /.tm line proved too sensitive to the def- 
inition of "pure" PAH 12.7 /im/11.3 /im ratio (e.g. as 
measured outside of the starburst cores where compar- 
atively little fine-structure line emission from H 11 re- 
gions is expected). Complications further arise from pos- 
sible intrinsic variations in the PAH ratios, extinction ef- 
fects (vli2.7^m/^ii.3/jin ~ 0.55), and unconstrained fine- 
structure line emission from disk Hii regions. The fluxes 
for [Ar ii] 6.99 ^m are much more reliable because possible 
contributions by other features in our sample galaxies are 
substantially smaller, as mentioned above. 



3.2. Mid-infrared images: mapping the features 

We obtained broad- and narrow-band images as well as 
maps of the PAHs and [Am] 6.99 /um line emission from 
the ISOCAM data cubes by applying to each pixel the 
procedures described above for the spectra. Figures ^, 
and present selected images and ratio maps for M82, 
NGC 253, and NGC 1808. Contours corresponding to the 
same levels relative to the peak intensity are plotted for 
all continuum and emission feature maps for ease of com- 
parison by visual inspection. A 3cr contour is also shown 
to delineate regions where observed small-scale structures 
are reliable. In the following, we describe the various maps 



source by source; their interpretation will be discussed in 
subsequent sections. 

3.2.1. M82 



The continuum and emission feature maps in M82 show a 
globally smooth spatial distribution, centered and peaking 
roughly 5" southwest of the nucleus. The PAH and 5.5 ^m 
continuum emission are the most extended and symmet- 
ric, with disk-like isophotes elongated along the galactic 
plane (P.A. « 70°). In contrast, the 15 /im continuum and 
[Am] 6.99 /im line emission have more compact distribu- 
tions which are more asymmetric relative to the nucleus. 
The 15/im/5.5/im ratio map outlines well the difference 
in peak morphology between the short and long wave- 
length continuum. The [Am] 6.99 /im distribution is the 
most compact, with centroid (determined from the emis- 
sion out to a radius of 25") displaced 3.5" southwest of 
that of the PAH and continuum emission, and showing 
only a slight extension towards the east. 

Our ISOCAM maps provide an important complemen- 
tary dataset to existing MIR images in the literature, 
which were mostly obtained in different bands or over lim- 
ited regions (although with higher angular resolution up to 
« 1"). The distributions observed in our PAH maps and 
for the PAH 3.29 /tm feature by Normand et al. (1995) 
and Satyapal et al. (1995) at « 1" resolution are consis- 



tent with each other. Our 15 /tm continuum map and the 
19.2 /im image of Telesco et al. (1991) are similar. Maps of 
the A^-band (10.8 /im) and 11.8 — 13.0 /tm emission gener- 
ated from the ISOCAM data cubes agree well with those 



of Telesco et al. ( |199lD and Telesco & Gezari ( |1992D within 
the regions covered by the latter two images. 
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Table 3. Intervals used to measure the continuum and features 
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Property 




Symbol 


Continuum points 
(/im) 


or intervals 


Intef 


^ration limits 
(/im) 


PAH emission " 




/pah 






6.00 


-9.00 


VSG emission 




/vsg 






13.5 


- 15.0 


5.5 /xm continuum 










5.40 


-5.52 


15.0 /im continuum 




./l5 f-ini 






14.8 


- 15.2 


PAH 6.2 fim 




F6.2 


5.81 - 5.99, 6.57 - 


6.80 


6.04 


-6.51 


PAH 7.7 /im 




F7.7 


7.14, 8.22-8.38 




7.19 


-8.17 


PAH 8.6 /im 




F8.6 


8.27 - 8.33, 8.84 - 


8.89 


8.33 


-8.84 


PAH 11.3 /tm 




i^ll.3 


10.95, 11.7-11.8 




11.1 


- 11.7 


PAH 12.7 /xm + [Neii] 


12.81 /im 


-F12.7 


12.1 - 12.2, 13.2 - 


13.3 


12.3 


- 13.2 


PAH 15.7 /tm + [Neiii 


15.56 fim 


i^l5.6 


15.0 - 15.2, 16.0 - 


16.1 


15.3 


- 15.9 


[Aril] 6.99 /im 




-PfArll] 


6.74 - 6.86, 7.14 - 


7.19 


6.86 


- 7.14 



" Integrated flux in the band including the short wavelength PAH emission complex. 

Integrated flux in the band probing the continuum emission from VSGs free from strong emission lines and PAH features. 



Table 4. Broad- and narrow-band measurements " 



Source 


Region 


/pah 


./vsg 




fl5 fim 


M82 


CAM FOV 


61.1 ±0.8 


78.3 ±0.3 


6.19 ±1.98 


80.8 ±0.6 




Disk 


12.7 ±0.1 


11.3 ± 0.1 


1.16 ±0.17 


10.8 ±0.1 




Core 


26.9 ±0.1 


45.7 ±0.3 


3.15 ±0.08 


49.3 ±0.5 




MIR Peak 


8.03 ±0.06 


16.3 ±0.2 


1.02 ± 0.04 


17.9 ±0.4 


NGC 253 


CAM FOV 


11.9 ±0.3 


25.3 ± 0.2 


2.12 ± 1.31 


27.2 ±0.3 




Disk 


2.49 ± 0.03 


3.93 ±0.03 


0.41 ±0.19 


4.12 ±0.06 




Core 


5.63 ± 0.03 


15.3 ±0.1 


1.05 ±0.11 


16.8 ±0.3 




MIR Peak 


1.68 ±0.02 


5.74 ±0.09 


0.35 ±0.02 


6.40 ±0.19 


NGC 1808 


Core 


2.66 ± 0.05 


3.13 ±0.04 




3.08 ±0.07 




MIR Peak 


1.10 ±0.01 


1.39 ± 0.03 


0.11 ±0.01 


1.41 ±0.05 



' AH flux densities are expressed in Jy. The uncertainties result from the formal effective uncertainties of the relative fluxes 
(see Sect. 



Table 5. Emission feature measurements ° 



Source Region 



Ft.- 



FSA 



Fi 



1.3 



F12. 



7 



Fl5 



F 



[Aril] 



M82 



CAM FOV 

Disk 

Core 

MIR Peak 



103 ± 1 
21.5 ± 0.2 

43.7 ±0.5 

12.8 ± 0.3 



226 ± 1 
46.5 ±0.2 
100 ± 1 
29.7 ±0.4 



31.8 ±2.4 
6.70 ±0.08 
13.8 ±0.4 
3.94 ±0.15 



57.7 ±0.9 
13.1 ±0.1 
22.0 ±0.5 
5.76 ±0.25 



54.0 ± 0.4 
9.74 ± 0.06 
26.4 ±0.5 
8.50 ±0.35 



5.77 ±0.47 
0.76 ±0.02 
3.49 ± 0.39 
1.43 ±0.21 



6.05 ±0.31 
0.67 ±0.02 
3.96 ±0.15 
1.64 ±0.08 



NGC 253 



CAM FOV 

Disk 

Core 

MIR Peak 



16.1 ± 1.1 
3.55 ±0.09 
7.30 ±0.11 
2.02 ± 0.06 



36.3 ±0.3 
7.63 ±0.07 
17.6 ±0.2 
4.95 ±0.09 



4.24 ±0.15 
1.03 ±0.02 
1.86 ±0.05 
0.54 ±0.02 



10.2 ±0.2 
2.41 ±0.03 
4.09 ±0.11 
1.07 ±0.06 



10.6 ±0.2 
2.02 ± 0.02 
5.71 ±0.17 
1.83 ±0.10 



0.55 ±0.12 
0.077 ± 0.007 
0.37 ±0.10 
0.19 ±0.05 



1.93 ±0.12 
0.36 ±0.01 
1.13 ±0.03 
0.37 ±0.02 



NGC 1808 



Core 

MIR Peak 



4.17 ±0.19 
1.59 ±0.05 



9.39 ±0.12 
3.74 ±0.08 



1.38 ±0.03 
0.57 ±0.02 



3.32 ±0.08 
1.31 ±0.05 



3.05 ± 0.06 
1.34 ±0.04 



0.096 ± 0.016 
0.039 ± 0.008 



0.23 ± 0.02 
0.087 ±0.007 



' All fluxes are expressed in 10 
fluxes (see Sect. 



Wm . The uncertainties result from the formal effective uncertainties of the relative 



Our [Am] 6.99 fim linemap globally resembles those of 
other tracers of ionized gas from H 11 regions at mid- and 
near-infrared wavelengths such as [Neii] 12.81 /im, [Arm] 
8.99 /im, [Siv] 10.51 /im, Bra, Br7, and Pa/3 (Achtermann 
Satyapal et al. 



& Lacy |I995 
1" 



1995). At a resolution of 



2", these maps reveal a rich sub-structure dominated 



by prominent sources « 6" and 12" southwest of the nu- 
cleus and « 6" to the nort heast (labeled Wl, W2, and 
El by Achtermann & Lacy 1995). While their intensity 
ratio depends somewhat on the emission line considered, 
Wl and W2 are together about three times brighter than 
El. Radial velocity data of the [Neii] 12.81 /im and Br7 
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emission (Larkin et al. 1994; Achtermann & Lacy 1995) 
suggest that most sources reside in a nearly edge-on ro- 
tating ring at radius coinciding in projection with Wl and 
El as well as along the stellar bar at larger radii, where 
the most recent starburst episode took place about 5 Myr 
ago (e.g. Forster Schreiber et al. ^002 ). Within the posi- 
tional uncertainties and resolution limitations, the spatial 
distribution of our [Aril] 6.99 /xm map peaks between Wl 
and W2 and encompasses El, and thus traces well the 
youngest starburst regions. 



Little differences are seen between the PAH maps 
but the PAH 6.2 ^m/7.7 iim and PAH 8.6 /zm/11.3 /xm ra- 
tio maps reveal spatial variations at the 20% and 60% 
level, respectively. The variations are significant within the 
brighter emission regions along the disk, where the relative 
uncertainties are < 15%. Structures seen towards the map 
edges are much less reliable as they become comparable in 
amplitude to the uncertainties. The PAH 6.2 /im/7.7 /im 
ratio is lower along the disk and reaches minima on each 
side of the nucleus. The overall morphology appears to 
curve northwards away from the nucleus and extensions 
are hinted at above and below the galactic plane. Similar 
shape and spurs are observed notably in the large-scale 
distribution of the molecular and ionized gas line emis- 
sion and of the radio continuum emission (e.g. Shen & 

1999| ). 



Lo |1995| ; Achtermann & Lacy |1995| ; Wills et al. 
The spatial variations in the PAH 8.6 /im/11.3 /im ratio 
match roughly those of the PAH 6.2 /im/7.7 /im ratio, with 
higher values along the disk and maxima flanking the nu- 
cleus. Noticeably, the western peak lies closer to the nu- 
cleus than the western PAH 6.2 /im/7.7 /im minimum, and 
the apparent curving and large-scale extensions have no 
counterpart in the PAH 8.6 /im/11.3 /im map. 



Figure ^ compares the PAH ratio maps wi th th e CO 
J : 1 ^ millimetric emission of Shen & Lo (1995) con- 



volved at the resolution of the ISOCAM maps. The over- 
all correspondence of the PAH 6.2 /im/7.7 /im minima and 
PAH 8.6 /im/11.3 /im maxima with the peaks in CO emis- 
sion, as well as the curved shape and northeastern exten- 
sion for the PAH 6.2 /im/7.7/im ratio, is quite striking. 
We believe that the observed variations in PAH ratios are 
mostly real. Characteristic patterns expected for ghosts 
are not seen (unresolved ring- or arc-like features most 
prominent in the presence of point-like sources). Artifacts 
due to the uncorrected flat field and straylight could pro- 
duce extended lobes on each side of an axis at P. A. w 150° 
for the M 82 data, i.e. roughly the minor axis (Biviano et 
al. [19984 |1998b| ; Okumura |2000D . However, such lobes 
would have a much larger extent than the features seen 
in our PAH ratio maps and the chromatic dependence be- 
tween 6.2 and 7.7 /im, and 8.6 and 11.3 /im is predicted 
and observed to be smaller than the measured variations 
at 20% and 60% levels, respectively. 



3.2.2. NGC 253 

In NGC 253, the emission in the continuum bands, PAH 
features, and [Am] 6.99 /im line is characterized by a 
strong peak within 1" — 2" of the nucleus, embedded 
in a diffuse envelope elongated along the galactic plane 
(P. A. Ks 50°). The 15 /im continuum and [Am] 6.99 /im 
line distributions do not seem to extend as far as the 
PAH emission in the outer parts of the source. The noisy 
A < 6 /im channels prevent reliable assessment of the 
lower level, large-scale 5.5 /im continuum emission. The 
centroids in the various images are essentially indistin- 
guishable (differences < 0.5 pixel). 

Though the images reveal little spatial structure be- 
cause of the intrinsically small source size and the lim- 
ited angular resolution, they are consistent with expecta- 
tions at this resolution based on previously published MIR 
maps (mostly obtained at ^ 1" resolution). These include 
broad-band images at 10.8 {N band) and 19.2 /tm, narrow- 
band images at 8.5, 10.0, 12.5, and 20.2 /im tracing PAH 
and/or continuum emission, maps of the 3.29 and 11.3 /im 
PAH emission, and images of the [Neil] 12.81 /im line and 
of the u nderlying con tinuu m + PAH emiss i on (P iiia et 
al. |l992l Telesco et al. Il993l Keto et al. |l993l |l999|: Kalas 



& Wynn- Williams |1994[ Boker et al. |1998|) We note that 
the [Nei l] linemaps of Boker et al. ( 1998|) and Keto et 
al. ( 1999 ) likely contain a contribution from PAH 12.7/im 
emission because the observations were made at low spec- 
tral resolution (FWHM w 0.2 /im). 

The spatial distribution seen in all our maps is ev- 
idently dominated by a very prominent compact MIR 
source, bette r out lined in arcsecond resolution images. 
Keto et al. ( 1999|) associated this source with a bright 
ofF-nucleus super star cluster resolved by optical Hubble 
Space Telescope observations. Based on their data, the 
source has a size of 20 pc (marginally resolved) and ac- 
counts for « 20% of the total continuum emission at 12 
and 20 /im, and 12% of the total [Neii] 12.81 /im flux. 
From our own measurements, nearly 25% of the total con- 
tinuum emission detected with ISOCAM between 13.5 
and 15 /im originates in the 70 pc-diameter MIR peak; 
this fraction varies between 10% and 30% for our various 
broad- and narrow-band and emission feature measure- 
ments (tables m and ||) . NGC 253 is quite remarkable in 
the compactness of its main MIR emitting region, with 
a FWHM ^ 150 pc while the optical disk extends over 
~ 20 kpc. Although of a different nature, the interact- 
ing system NGC 4038/4039 offers a similar example, with 
15% of its 12.5 — 18 /im luminosity being produced within 
a 100 pc-size star-forming knot in t he ove rlapping region 
between the galaxies (Mirabel et al. 1998). 

Additional small-scale structure of NGC 253 at MIR 
wavelengths includes a secondary much fainter peak in 
the 12 and 20 /im continuum emission nearly coinciding 
with the nucleus, about 2" northeast of the prominent 
source discussed above. The [Neii] 12.81 /im line emission 
differs somewhat from the continuum, except possibly for 
the brightest peak, showing a bilobal or arc-like structure 
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also hinted at in Br7 images and suggestive of a circum- 
nuclear star-forming ring (Boker et al. 199S; Engelbracht 
et al. 1998; Keto et al. 1999). None of these features, how- 
ever, is resolved with ISOCAM. 

The PAH 6.2 /im/7.7 fim and PAH 8.6 ^m/11.3 fim ra- 
tio maps show variations of about 40% and nearly a 
factor of 2, respectively. These are statistically signifi- 
cant in view of the corresponding formal uncertainties 
of < 20% and < 30%. The ratio maps differ markedly, 
with the PAH 6.2 /im/7.7/im image indicating a general 
increase from south to north of the nucleus whereas the 
PAH 8.6 /im/11.3 /im ratio appears more centrally con- 
centrated, peaking near the nucleus and slightly more ex- 
tended northeast. The reality of the arc-like feature « 10" 
northwest of the nucleus in the PAH 8.6 /xm/11.3 /xm map 
is dubious because it is barely resolved and could be due 
to ghost effects given the strong unresolved peak of the 
emission. No corresponding structure that could perhaps 
support this feature is seen in maps of the molecular gas 
emission, of the radio continuum emission, of tracers of 
ionized gas or even of extinction (e.g. Peng et al 



1996 



Ulvestad & Antonucci |1997| ; Engelbracht et al. |1998| ). The 
limited region for which S/N > 3 prevents us from exam- 
ining the larger-scale distribution in PAH ratios. 
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3.2.3. NGC 1808 

The starburst in NGC 1808 covers a region of comparable 
physical size to that of M 82 but being three times more 
distant, less structural details are resolved by ISOCAM. 
The PAH emission appears the most extended and ori- 
ented parallel to the major axis of the galaxy (P. A. « 
140°), peaking at the nucleus. The 15 ^m emission fol- 
lows closely the PAH emission. The [Ar ii] 6.99 fj,m line 
emission is the most distinct in that the peak clearly is off- 
nucleus, about 5" to the southeast, and the centroid of the 
emission region lies 2.5" southeast of that for the 15 /^m 
continuum and PAH emission. The [Am] 6.99 /im map 
agrees well in peak position and extent with the global 
distribution of the most intense star-forming regions, or 
hot spots, as traced by H recombinatio n line s and radio 
continuum emission (e.g. Saikia et al. 199C; Krabbe et 
Kotilainen et al. 



al. 1994 



19961 ). Due to the limited num- 



ber of pixels with S/N > 3, no useful ratio maps could be 
made for NGC 1808. 



Previously published 
(Telesco et al 



1993 



N band image s of NGC 1808 
Krabbe et al. 2001) show an over- 



all similar morphology as our 15 /im continuum and PAH 
feature maps with, at « 1" resolution, a strong point- 
like source at the nucleus and a southeastern exten- 
sion covering the starburst regions. ISOCAM broad-band 
LW4 (6 /ini) polarisation observations were presented by 
Siebenmorgen et al. (2001) along with a CVF spectrum 
of the central 25" of NGC 1808 which is essentially iden- 
tical in shape to ours of the starburst core (30" aperture; 
Fig. ^). Siebenmorgen et al. (2001) assigned all detected 
features to PAHs and successfully reproduced them with 



Fig. 8. Comparison of ISOCAM PAH ratio maps of M 82 
with the distribution of CO J : 1 ^ millimetric emis- 
sion. The PAH 6.2^m/7.7^m and PAH 8.6 /xm/11.3 /xm 
maps (top and bottom panels, respectively) are displayed 
as grayscale images as for Fig. |[ The contours sho w the 
CO J : 1^0 emission from the map of Shen & Lo ( 1995[ ) 
convolved at the resolution of the ISOCAM images (from 
the original 2.5" to 5.6"). The physical linear scale and the 
FWHM of the ISOCAM PSF are indicated at the bottom 
left and right. The cross at relative coordinates (0", 0") 
marks the location of the galaxy nucleus. 



PAH emission alone in their radiative transfer models. As 
also noted by these authors, the contribution of nebular 
gas emission lines blended with PAH features is uncer- 
tain for NGC 1808, and we cannot constrain it using, e.g., 
SWS data as for M82 and NGC 253. However, our [Am] 
6.99 map clearly differs from those of the PAH emis- 
sion and is consistent with the spatial distribution of the 
bulk of Hii regions, supporting the idea that the 7 fim 
feature is indeed dominated by the [Am] line instead of 
PAHs. 



3.3. Salient features 

We summarize here the most important aspects of the 
results presented above: 

— in terms of features present, the MIR spectral energy 
distribution is essentially identical among the three 
galaxies as well as within each of them. 
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— the largest spectral variations are observed in the 
A ^ 11 /im continuum intensity relative to the shorter 
wavelength emission, while the shape of the A = 5 — 
11 /im region is nearly invariant, 

~ the PAH emission is spatially more extended than the 
15 /im continuum and [Aril] 6.99 /im line emission, 

— the [Am] 6.99 /xm line emission traces well the young 
starburst regions, and 

— the relative intensities of the PAH features can exhibit 
complex spatial variations. 

4. Origin and spatial distribution of the PAH and 
continuum emission 

The ubiquity of PAH features - seen towards sources 
as diverse as circumstellar regions, diffuse atomic clouds, 
Hll regions, molecular clouds, normal galaxies, starbursts, 
and ULIRGs - and the close association of the MIR 
VSG continuum emission with Galactic and extragalac- 
tic actively star-forming sites now seem well established; 
the ISO mission played a major role in this recognition 
(see Geballe |1997[ Tokun aga |l997t Tielens et al |1999 ; 
Cesarsky & Sauvage 1999 , Genzel & Cesarsky 2000 for re- 
views). Spatial mapping of Galactic sources has revealed 
in more detail the origin of the main MIR emission com- 
ponents in star-forming regions: PAH features arise pre- 
dominantly in the photodissociation regions (PDRs) at 
the interface between H ii regions and molecular clouds 
while the VSG A <; 11 /im continuum is more intense and 
steeper closer to the exciting source (e.g. Verstraete et 
al. |1996| ; Cesarsky et al. |19964 |1996h| ; Crete et al. |1999D . 
By extension, in star-forming galaxies, the emission com- 
ponents are believed to trace PDRs and Hii regions, re- 
spectively. This is suppor ted by spectral deco mposit ion of 
starburst galaxies (Tran 1998 ; Laurent et al. 2000 ), and 
by photometric and spcctrophotometric imaging of spiral 
galaxies and starburst systems (e.g. Mirabel et al. 1998 ; 
Roussel et al. 2001a ). Low excitation fine-structure line 
emission such as [Am] 6.99 /im, [Ami] 8.99 /im, [Neii] 
12.81 /im, and [Neiii] 15.56 /im usually accompanies the 
long- wavelength VSG continuu m and originates primarily 
in Hii regions (e.g. Sturm et al. 200C). 

The interpretation of the ISOCAM maps of M82, 
NGC 253, and NGC 1808 within the above framework is 
not straightforward. Figure ^ complements the various 
maps with light profiles of the 6 — 9 /im, 13.5 — 15 /im, and 
[Am] 6.99 /im emission taken along the major and minor 
axes of the galaxies. The profiles clearly support the exten- 
sion along the major axis indicated by the maps, and show 
that the emission in these tracers along the minor axis is 
but marginally resolved and essentially indistinguishable 
within each galaxy. 

Somewhat surprisingly given their common origin in 
Hii regions, a close spatial correlation between the VSG 
continuum and the [Am] 6.99 /im line emission is not 
observed in our sample galaxies. The exception may be 
NGC 253 but this probably results from the intrinsi- 
cally small MIR source size and resolution limitations 



(Sect. 3.2.2| ). Although the extension to the northeast is 
real, the differences between the various profiles even along 
the major axis are smaller than the PSF FWHM. We 
conclude for NGC 253 that on scales comparable to and 
larger than the resolution of our ISOCAM maps (5.2", or 
~ 60 pc), no significant differences are observed in the 
spatial distribution of the PAH, VSG, and [Am] 6.99 /im 
emission within the central 400 pc. 

For M82, the [Am] 6.99 /im emission along the galac- 
tic plane peaks at a different position and is more compact 
than the VSG and PAH emission. The different distribu- 
tions of the VSG and [Am] 6.99 /im emission could be 
interpreted in terms of differences in average energy of 
the exciting photons: the [Am] 6.99 /im line traces the 
> 16 eV radiation field while the VSGs may be heat ed by 



UV photons at lower en ergie s as well (Desert et al . [1990 



Siebenmorgen & Kriigel 1992; Boula nger et al. 1994; Dwek 
et al. [I997t Jones & d'Hendecourt |2000| ). The brighter 
PAH emission around the nucleus could be originating 
predominantly in the PDRs associated with the Hil re- 
gions within the most active starburst sites. Beyond a ra- 
dius of « 300 pc, the PAH profile is resolved out from the 
VSG profile and more extended. These differences could 
be attributed to excitation of aromatic band carriers by 
yet softer photons or to an increased filling factor for the 
PDRs compared to the H II regions in less intense starburst 
regions, and perhaps to a contribution from diffuse cirrus 
clouds as observed in the Milky Way and in some spiral 



galaxies (e.g. Ristorcelli et al. 1994; Onaka et al. 
Mattila et al. 



1999; Roussel et al. 2001a; Li & Draine 



1996 



2002). Alternative possibilities could include destruction 
by photodissociation or depletion of the ISM in PAHs and 
VSGs in the most intense starburst regions, flattening the 
sp atial d istribution around the peak (e.g. Boul anger et 
al. |1988| ; Carral et al. |l994| ; Normand et al. |l995|) . 

In NGC 1808, the bulk of the [Am] 6.99 /im emis- 
sion is clearly displaced to the southeast compared to 
the emission in the PAH and VSG bands defined above, 
which h ave v irtually identical distributions. As discussed 
in Sect. 3.2.5 , the [Am] 6.99 /im distribution agrees well 
with the off-nucleus location of the main star-forming 
sites. The spectrum integrated over the starburst core and 
MIR peak in NGC 1808 is also much flatter than observed 
towards the corresponding regions in M 82 and NGC 253 
and resembles more the spectrum of their outer disks 
(Fig. |). The long- wavelength continuum in NGC 1808 
may in fact not be dominated by VSG emission from 
pure H II regions but rather produced primarily in PDRs 
by small particles akin to the carriers of the main UIB 
bands. This is reminiscent of the situation in disks of spi- 
ral galaxies. There, the 15 /im emission correlates with 
the 7 /im PAH-dominated emission and lies in a distinct 
regime compared to more active regions including circum- 
nuclear regions in spiral (especially barred) galaxies and 
starbursts, characterized by an excess in the 15/im/7/im 
colour (e.g. Roussel et al. 2001a; Dale et al. 2001; Forster 
Schreiber, Roussel, & Sauvage, in prep.). From radiative 
transfer modeling of NGC 1808's ~ 3 — 200 /im spectrum, 
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Siebenmorgen et al. ( |200l| ) concluded that the MIR range 
is dominated by PAHs and large dust grains (radii be- 
tween 100 and 2400 A) with negligible contribution by 
very small graphites and silicates (radii of 5 — 80 A). The 
key feature of their models is the inclusion of hot spots 
where large grains are heated locally by massive stars to 
higher temperatures giving rise to the 25 /im emission and 
gradually contributing less towards shorter wavelengths 
where PAHs take over as main emitters. 

We note that the morphology of the [Am] 6.99 /xm 
emission for M 82 and NGC 1808, especially the off-nucleus 
peak position, reflects the spatial distribution of the H ii 
regions and not variations in the excitation state of 
the ionized gas whereby, e.g., [Arm] emission would be- 
come important at the expense of [Ar ii] . As discussed in 
Sect. 3.2.1 , the Hii regions in M82 are concentrated in 
a circumnuclear ring- like structure at radius of « 85 pc. 
Direct tracers such as near-infrared to radio H recombi- 
nation lines, MIR Ne, Ar, and S fine-structure lines, the 
He I 2.06 fj,m recombination line, and radio thermal contin- 
uum all show a similar distribution characterized by peaks 
flanking the nucleus, with the southwestern complex being 
the most intense, and little emission at the nucleus (e.g. 
Satyapal et al. 1995| ; Achtermann & Lacy 1995| ; Seaquist 
et al. [1996 ; Forster Schreiber et al. 2001). In particular, 
the [Ann] 8.99 /im map of Achtermann & Lacy (1995) 
tracing the > 28 eV radiation field clearly shows this 
bilobal distribution and relative lack of emission at the 
nucleus. Less information is available for NGC 1808 but 
Br7 imaging clearly reveals the H ii region complexes to 
lie predominantly to the east and southeast of the nucleus, 
thus indicating a circumnuclear location for the starburst 



1996) 



(Krabbe et al. 1994; Kotilainen et al. 

Extinction effects could partly account for differences 
in the PAH, VSG, and [Am] 6.99 /im emission on small 
spatial scales (comparable to the resolution elements in 
our ISOCAM data) but are not likely to affect the rela- 
tive concentrations. The highest levels of obscuration are 
measured in the near surroundings of the MIR peaks in all 
three galaxies and gene rally decrease at la rger d istances 
(e.g. Larkin et al. |1994| ; Engelbracht et al. |1998| ; Krabbe 
et al. 1994). Depending on the extinction law applicable 
in the 3 — 10 nm region (Sect. 5.1 below), such a radial 
gradient in extinction would cither affect very little the 
relative global distributions of the tracers discussed here, 
or suppress more severely the shorter wavelength emis- 
sion in the central regions, which is inconsistent with the 
overall trends for the 15 /im continuum and [Ar ii] 6.99 /im 
line. 

The nature of the short-wavelength continuum in nor- 
mal and pure starburst galaxies is still debated. It is much 
weaker than in systems hosting an active galactic nucleus 
(AGN) where it is attributed to hot dust ( 150 - 1700 K) 
close to the AGN (Genzel & Cesarsky 2000 and references 
therein). In non-AGN galaxies, the existence of 3 — 5 /im 
featureless continuum from small transiently heated dust 
grains underlying the PAH features has been proposed 



tures are best fitted with Lorentzian profiles whose wide 
wings can account for the apparent continuum pedestal 
(e.g. Boulanger et al. 1998b). Our narrow-band 5.5 /xm 
continuum map of M 82 provides only tentative evidence 
of a more extended distribution than for the 15 /im con- 
tinuum emission which would suggest a closer association 
of the emitting particles with PAHs than VSGs (see also 
Mattila et al. 



1999) 



(e.g. Helou et al. 2000). On the other hand, the PAH fea- 



5. PAH and continuum emission: relations to the 
ISIVI conditions and star-forming activity 

The two most striking aspects seen in the spectra pre- 
sented in Fig. 1^ are their remarkable similarity in the 
A = 5 — 11 /xm regime dominated by PAH emission, and 
the spread in relative intensity of the longer wavelength 
emission attributed to VSGs. In this section, we focus on 
these two spectral components after discussing extinction 
effects at MIR wavelengths. 

5.1. The effects of extinction at mid-infrared 
wavelengths 

Interstellar extinction can significantly affect the MIR 
spectral energy distribution (SED) in obscured sources. 
In this regard, Rigopoulou et al. (1999) noted the grad- 
ual suppression of the PAH 8.6 /im feature in progressively 
obscured starbursts and ULIRGs. They further proposed 
that variations in the PAH 6.2 /im/7.7 /im ratio are dom- 
inated by extinction effects on the basis of the trend of 
decreasing ratio with increasing extinction observed in a 
subset of their sample for which reliable extinction deter- 
minations were available. 

Another important consideration concerns the trough 
centered near 9.7 /im traditionally attributed to absorp- 
tion by interstellar silicate dust grains. Recent studies have 
cast doubt on the reliability of extinction estim ates b ased 
on the observed feature depth. Sturm et al. (2000) dis- 
cussed this issue in detail based on SWS data of M82 
and NGC 253, from a comparison with the spectrum of 
the Galactic reflection nebula NGC 7023 and from con- 
siderations of the relative optical depths expected for sil- 
icate absorption near 9.7 and 18 /im. In particular, they 
showed that the M 82 SWS spectrum, including the dip at 
9.7 /im, can be well reproduced by a combination of the 
NGC 7023 spectrum and a power law rising from 8.5 /im 
longwards without invoking large extinctions. More gen- 
erally, the overall spectral invariance in the 5 — 11 //m re- 
gion among a variety of Galactic and extragalactic sources 
which are known or expected to cover a large range in ex- 
tinction seems to argue in favour of the 9.7 /im "absorption 
feature" being predominantly due to the gap between the 
main PAH emission complexes (e.g. Helou et al. 2000). 

These suggestions raise the important question of how 
such an interpretation can be reconciled with the large 
and variable extinction in M 82, NGC 253, and NGC 1808 
determined by alternative methods based, for instance, 
on the relative intensities of H recombination lines. To 
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address this issue, we explored quantitatively the effects 
of extinction at MIR wavelengths by applying a range 
of extinction to a representative template spectrum con- 
structed in a similar way as Sturm et al. (2000). We 
added to a PDR component (from the ISOCAM data of 
NGC 7023) a power law fi, (x {X — 8.5)" which represents 
well the VSG emission seen in our data. We adopted an in- 
dex a = 1.5 and scaled the two components so as to obtain 
resulting spectra similar to that of M82. For the purpose 
of illustrating the effects of different levels of extinction, 
the choice of these parameters is irrelevant. 

We considered two representative geometries: a uni- 
form foreground dust screen and a homogeneous mixture 
of dust and sources. For these models, the observed in- 
tensity /(A) of the emerging radiation is proportional to 
the total intrinsic intensity of the sources Iq (A) by the fac- 



tors e "^^ and (1 



^ ) /tx , respectively, where t\ is the 



optical depth related to the extinction in magnitude via 
A\ = 1.08 6 TA. We adopted the extinction law proposed 
by Draine ( 1989| ) and investigated the effects of deviations 
from this law in the 3 — 10 /im region as found towards the 
Galactic Center (Lutz 1999; hereafter simply "GC law"). 
Such deviations are consistent with the H recombination 
line spectrum observed with the SWS in M82 (Forster 
Schreiber et al. 2001). We varied the level of extinction 
specified in visual magnitudes Ay for the mixed model in 
the range Ay^-^ = — 700 mag. For a meaningful com- 
parison, we varied Ay for uniform foreground obscuration 
in the range 

^UFS ^0-50 mag, which results in the 
same attenuation at 10 fim as the range considered for the 
mixed model. We normalized each extincted model spec- 
trum to the integrated flux between 6.0 and 6.6 /im, as for 
the spectra of Fig. |[ 

Figure |l^ presents our simulations. The plots outline 
the distinct behaviours for the two geometries considered 
with increasing Ay. They especially emphasize the fact 
that for the mixed model, differential extinction effects 
over the spectrum rapidly reach an asymptotic regime 
where large variations in Ay become less and less percep- 
tible as Ay increases. In addition, except for the most ex- 
treme cases {Ay > 20 mag) , the choice of dust and source 
geometry has little impact. This is because the absolute 
levels of extinction are fairly small and the wavelength de- 
pendence of extinction is relatively weak; the wavelength 
range covered is thus too limited to probe appreciably dif- 
ferent optical depths which could allow the discrimination 
between different dust and sources distributions. 

For the Draine extinction law, the effects are largest 
between 8 and 13 /im. Our models illustrate well how the 
PAH 8.6 and 11.3 /im features are substantially suppressed 
as extinction increases, indicating that their ratio with the 
PAH features at 6.2 and 7.7 /im can be much affected by 
obscuration. Moreover, for a continuum level defined over 
the limited 8 — 13 /^m interval, relatively small optical 
depths would be inferred even if the extinction is in fact 
large. Assuming the GC law at 3 — 10 /.im results in a 
remarkably different behaviour. The extinction effects are 
significantly smaller in the 8 — 13 /xm region and the sup- 



pression of the PAH 8.6 and 11.3 /xm features is much less 
important. The largest effects are observed in the relative 
level between the short and long wavelength regions. This 
is a consequence of the much flatter GC law between 3 
and 10 /xm. 

Our simulations indicate that dust obscuration adds 
a significant degree of degeneracy in the interpretation of 
the observed MIR emission, complicating the determina- 
tion of intrinsic properties. Depending on the extinction 
law used, the shape of the PAH 8.6 and 11.3 /xm features as 
well as their flux ratio with the PAH 6.2 and 7.7 /xm fea- 
tures vary considerably. Therefore, extinction estimates 
based on diagnostics involving the PAH 8.6 or 11.3 /xm 
feature may be strongly biased by the choice of extinction 
law. Furthermore, part of the spread in intensity observed 
in the long-wavelength continuum relative to the shorter 
wavelength emission could be attributed to varying ex- 
tinction levels (see Fig. I and Sect. U below). 

The interpretation of the 9.7 /xm trough as largely due 
to the gap between the flanking PAH complexes appears 
to hold over a wide range of extinction, especially for the 
GC extinction law. The near-invariance of the 5 — 11 /xm 
spectrum is therefore not inconsistent with large and/or 
variable obscuration among and within galaxies. We il- 
lustrate this with the case of M82, for which a mixed 
model with Ay = 52 mag best reproduces the observed 
H recombination lines from radio to optical wavelengths 
while the best-fit uniform foreground screen extinction of 
Ay ~ 5 mag provides a much poorer fit (Forster Schreiber 
et al. ^001 ). We applied these two extinction models to the 
same template as for the simulations described above. The 
resulting spectra are plotted along with the observed spec- 
trum of M 82 in the bottom panels of Fig. |l^. We did not 
attempt to fine-tune the models by formal fitting. Given 
the uncertainties on the exact nature of the emitting par- 
ticles, specific assumptions on model parameters are not 
well constrained and a simple empirical approach is suffi- 
cient for our purposes. The comparison shows that strictly 
from the point of view of the 5 — 16 /xm range, both ex- 
tinction models reproduce equally well the observed SED 
of M 82, thus demonstrating that a high extinction cannot 
be excluded from the overall shape of the MIR spectrum 
alone. 

We wish to emphasize that we do not dismiss the pos- 
sibility of silicate absorption around 9.7 /xm in general but 
merely want to point out the difficulties involved in the 
interpretation of the observed feature. The SEDs of mod- 
erately to highly obscured sources including ULIRGs of- 
ten exhibit a strong dip near 9.7 /im together with sup- 
pressed PAH 8.6 and 11.3 /im features relative to those 
at 6.2 and 7.7 /im, consistent with the presence of silicate 



grains (e.g. Dudley & Wynn-WiUia ms |1997| ; Dudley |1999 ; 
Laurent et al. 2000; Le Floc'h et al. ^002 ). Given the near 



invariance of the PAH spectrum in a wide range of en- 
vironments, it should be possible to define an indicator 
measuring the differential extinction between the 9.7 /im 
region and adjacent less affected intervals that quantifies 
the absolute extinction. One must however be aware of the 
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importance of sufficient wavelength coverage and resolu- 
tion to assess properly the impact of the PAH complexes 
which can make the silicate absorption look artificially 
deep, of the dependence of the inferred Ay on the assumed 
extinction law and spectral intervals used to measure the 
feature depth, and of the limited sensitivity of this diag- 
nostic leading to rather large uncertainties in the derived 
Ay - The upcoming launch of SIRTF in 2003 will provide 
5 — 40 ^m spectroscopy with an increase in sensitivity by 
two orders of magnitude and will help us address this issue 
by enabling a better sampling of the continuum emission 
as well as both the 9.7 and 18 /im silicate bands. 



5.2. The \> 11 fim continuum emission 

The A ^ 11 /im continuum emission in our data exhibits 
a large spread in intensity relative to the shorter wave- 
length emission. Similar variations have been seen in nor- 
mal spirals and starburst galaxies observed with ISO in- 
struments and, based on evidence provided by Galactic 
sources, are generally interpreted in terms of the rela- 
tive contribut ion of H ii regions t o the M IR emissio n (e.g . 
Laurent et al. |2000| ; Roussel et al. |200l4 Dale et al. |200lD . 
A close link between the 15/xm/7/im colour as measured 
through the ISOCAM broad-band LW3 and LW2 filters 
(A = 12 — 15 /im and 5 — 8.5 /im, respectively) and the 
IRAS 25 /tm /12 /im colour has been em phasized by Dale 
et al. ( |200lD and Roussel et al. ( |2001a{ ). 

In order to assess quantitatively the relationship be- 
tween the long-wavelength continuum properties and the 
H II regions, we investigated the relationship between the 
15 fim narrow-band continuum and the [Ar ii] 6.99 /im line 
emission. The [Am] 6.99 /im is the most direct probe of 
H II regions and the least contaminated by PAH emission 
available from our ISOCAM data sets. The Ar" ionization 
potential of 15.8 eV is close to that of H (13.6 eV). Argon 
being a noble element, it is, like neon, not expected to 
be significantly depleted onto dust grains; while the argon 
gas-phase abundance may increase with time as a result of 
star formation activity and differ among galaxies, it is not 
likely to vary significantly over scales of ~ 100 — 1000 pc 
as covered by the starburst regions in our sample. For 
the relatively low nebular excitation in all three sources 
(Fi g. 1^; see also Sturm et al. 2000 ; Forster Schreiber et 
al. ^OOlD , argon is mostly singly-ionized so that the [Ar ii] 
6.99 /im emission should trace the Hii regions in much 
the same way as H recombination lines and its luminosity, 
scale with the radiation field intensity of the young stars. 
In addition, since the gas-phase abundances of the galax- 
ies are solar within a factor of a few (e.g. Webster & Smith 
1983| ; Forbes et al. |1992| Forster Schreiber et al. ^OOlD , the 



proportionality factor should be roughly similar. 

Figure ^ presents the 15 /im continuum versus [Ar ii] 
6.99 fim luminosities for our galaxies, normalized to unit 
projected surface area (Sis^m and Ejai II])- The data are 
shown for the selected regions as well as for the individ- 
ual resolution elements enclosed within the radii defining 



the outer limit of the disk regions in M82 and NGC 253. 
Resolution elements with measurements at < 3(7 were ex- 
cluded. The data points follow a well defined distribution, 
with a remarkable overlap for the different galaxies. Least- 
squares fitting to the resolution elements' data accounting 
for the individual formal uncertainties yields the relation 



l0g(I]i5^m) = 1.01 l0g(S[A,. II] ) +0.807, 



(1) 



with dispersion around the fit of 0.203 dex. The extinc- 
tion paths in the diagram of Fig. ^ are nearly parallel to 
the locus formed by the data points; however, implausibly 
large differential extinctions would be required in order to 
account entirely for the observed distribution. 

The essentially linear proportionality we find between 
Sis^m and S[Arii] indicates that the 15 /im emission pro- 
vides a good quantitative indicator of the star-forming 
activity in starburst environments, to within uncertainties 
determined by the dispersion of the data and by extinction 
(a factor of 2.5 for an Ay = 50 mag assuming purely fore- 
ground extinction). Roussel et al. ( 2001b ) reached a sim- 
ilar conclusion for more quiescent spiral disks where the 
12 — 15 /im emission scales linearly with the Ha line flux. 
Part of the scatter in Fig. |ll] may result from variations 
in the physical conditions and exact composition of the 
gas and dust within and between the galaxies. More im- 
portant factors, however, are variations in relative spatial 
distribution of the emission and possibly of the excitation 
state of the gas. 

The observed relationship holds remarkably well in 
view of the differences in morphology of the 15 /im contin- 
uum and [Am] 6.99 /xm line emission as seen in Figs. |], 
^, 1^, and ||. The correlation extends however over two to 
three orders of magnitude, much larger than the disper- 
sion of a factor of 1.6. Undoubtedly, the relative spatial 
variations between both tracers contribute significantly 
to the scatter. Within the starburst cores of M82 and 
NGC 253, where the steeply rising SEDs at A ^ 11 /im 
are dominated by VSG emission, the differences may be 
attributed to different ranges in exciting photon energies 
for VSGs and [Am] 6.99 /xm line as well as to extinc- 
tion effects (Sect. |). For NGC 1808, the spatial distri- 
butions differ probably because a different dust/particles 
population produces the flat long-wavelength continuum 
(Sect. 1). Nevertheless, the data for NGC 1808 and for the 
disk regions of M82 and NGC 253, characterized by flat 
^11 /im SEDs, are well described by Eq. |l|. 

Variations of the excitation state of the photoionized 
nebulae may influence the Sis^^m versus S[Ai II] relation- 
ship as well. In particular, the [Ami] 8.99 /im line is fairly 
strong in the ISO-SWS spectrum of M82 (Fig. |) but 
the [Arm] 8.99/im/[Am] 6.99 /im ratio of 0.18 (0.26 af- 
ter extinction correction) is low and abundance estimates 
indicate that ^ 25% only of the argon is doubly ion- 
ized (Forster Schreiber et al. pOOl ). As emphasized in 
Sect. ^ the spatial distributions of various MIR fine- 
structure lines and Bra are similar and suggest a roughly 
constant excitation state of the Hii regions across M82 
(Achtermann & Lacy 1995). This is further confirmed by 
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Fig. 11. Relationship between MIR tracers of H II regions 
in M 82, NGC 253, and NGC 1808. The diagram shows the 
variations of the 15 /xm continuum luminosity as a func- 
tion of the [Aril] 6.99 fim line luminosity normalized per 
unit projected area. Data are plotted for selected regions 
(open symbols labeled "F" for the ISOCAM field of view, 
"D" for the disk regions, "C" for the starburst core, and 
"P" for the MIR peak). The filled circles and crosses show 
the values for individual resolution elements for M 82 and 
NGC 253, respectively, within radii of 30" and 15" (outer 
disk annulus; tabled); measurements at < 3(t arc excluded 
(4 points for M82, II for NGC 253). Typical formal un- 
certainties are indicated by the error bar at the bottom 
left and are smaller than the symbol sizes; they are 5% 
for Eis^m and 10% for S[Ai II] (the median and average 
differ by < 1.5%). The straight line is a least-squares fit to 
the resolution elements' data accounting for the individ- 
ual uncertainties. The effects of extinction are shown in 
the right part of the plot for a uniform foreground screen 
and a mixed model ("UF S" and "MIX"), for the extinc- 
tion law of Drainc (1989, "D89") and with modifications 
at 3 — 10 /im as found towards the Galactic Center (Lutz 
I99£, "GC"). Selected values of extinction in visual mag- 
nitudes Ay are labeled. 



the nearly uniform Hei 2.06/im/Br7 ratio in the central 
260 X 160 pc which corresponds closely to the excitation 
derived from the SWS [Neiii] 15.56 ^m/[Ne 11] 12.81 ^m, 
[Arm] 8.99^m/[Arii] 6.99 /^m, and [Siv] 10.5 ^m/[S ill 
I8.7/im ratios within larger apertures up to 430 x 225 pc 
(Forster Schreiber et al. ^001 ). Similar arguments are 
more difiicult for NGC 253 and NGC 1808 because of the 
lack of relevant data. However, the [Arm] 8.99//m/[Arii] 
6.99 fiui Ri 0.03 measured from the SWS spectrum of 
NGC 253 is substantially lower than for M82, as is the 
[Neiil] 15.56 ^m/[Ne 11 ] 12.8 1 ^m ratio (0.06 compar ed to 
0.18; Thornley et al. pOOO| see also Giveon et al. p002 



for the relationship between these two MIR line ratios). 
For NGC 1808, we can only note from the ISOCAM spec- 
tra (Fig. H) that the weakness of the 15.7 /xm feature 
and the non-detection of [Arm] 8.99 /Ltm indicate low ex- 
citation supporting that argon mostly is singly ionized. 
Summarizing, the excitation state of the gas might cause 
for M 82 a general but small offset to the left in the rela- 
tionship of Fig. |ri| compared to the other galaxies because 
of the larger fraction of Ar while spatial variations are 
not likely to introduce scatter larger than this offset. These 
effects are expected to be smaller for both NGC 253 and 
NGC 1808. 

The relationship between the 15 /im continuum and 
star formation intensity is supported by the analysis of 
a larger sample including spiral and starburst galaxies 
that will be presented in a forthcoming paper (Forster 
Schreiber, Roussel, & Sauvage, in prep.). In this paper, 
based on the previous work of Roussel et al. ( |200Ib ) for 
spiral disks, we find that the 15 /im continuum as well as 
the 7 /im PAH-dominated emission correlate well with the 
star formation rate over nearly 7 orders of magnitude, with 
the interesting distinction of two regimes that correspond 
to quiescent star formation in disks and more intense ac- 
tivity in circumnuclear regions and starbursts. 



5.3. The X = 5 — 11 emission 

The near invariance of the 5 — II fim spectrum seen in our 
data of M82, NGC 253, and NGC 1808 has been noted in 
a number of other studies of different types of galaxies 
powered mainly by star formation, as well as in a variety 
of Galactic sources (see Tielens et al. 1999 for a review 
see also e.g. Boulanger et al 



|l998a| ; Helou et al. pOOO ; 
Uchida et al. 2000[ ). This indicates that the PAHs are very 
stable under a wide range of physical conditions despite 
their small sizes (typically ~ 100 atoms) and is taken as 
direct observational evidence for the stochastic nature of 
the emission processes involved. Furthermore, if different 
types of PAHs coexist, the near constancy in the 5 — 11 /zm 
region suggests that their relative abundances vary little. 
Finally, it has consequences on the interpretation of and 
extinction measurements from the 9.7 /im silicate absorp- 
tion feature because of the intrinsic gap between the main 
6 — 9 /im and II — 13 /im PAH complexes (Sect. 5.1). 

Variations of the relative intensities of the 
PAH features do exist, however. We focus on the 
PAH 6.2 /im/7.7 /im and PAH 8.6 /im/II.3 /im ratios, 
maps of which were presented for M 82 and NGC 253 in 



] Sect. 3.2. Of particular interest is also the enhancement 



of the PAH II.3/im relative to the other features in 
NGC 1808 compared to M82 and NGC 253, evident 
in Fig. 1^. The variations of PAH ratios in our sample 
are comparable to those observed within and between 
Galactic sources of similar types (e.g. Cesarsky et 
al. |I996at Boulanger et al. |I998a|; Lu |I998| Crete et 



al. 1999; Uchida et al. 2000), along the disk of the spiral 
galaxy NGC 891 (Mattila et al. |I999D , and among a 
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sa mple of 15 starbursts and ULIRGs (Rigopoulou et 
al. 1999 ). A detailed discussion of PAH ratios is beyond 
the scope of this paper, but we briefly mention possible 
interpretations of the variations seen in our data in the 
light of some recent theoretical and empirical work. 

Since the emission mechanism is stochastic in na- 
ture, PAH ratios are not directly related to the SED 
of the incident radiation (e.g. Boulanger et al. 1998a; 
Uchida et al. 200C). On the other hand, they may 



carry indirect information since PAHs exposed to hard 
and intense radiation fields can be ionized, lose H 
atoms, or be photodissociated (e.g. Leger et al 



198S 



Allam andol a et al. 1989| ; Schutte et al. 1993|; AUamandola 



et al. |199£ ; Hudgins & Allamandola 1999). In partic- 
ular, the PAH 8.6 /im/11.3 ^m is believed to trace the 
fraction of singly ionized to neutral PAHs, presumably 
driven by the strength of the radiation field fro m th e 
OB stars (see refere nces above, and Joblin et al. 1996| ). 
Draine & Li ( 2001 ) presented a thorough calculation 
of the expected PAH spectrum as a function of vari- 
ous parameters including PAH size, charging conditions, 
and starlight intensity. Their results indicate that the 
PAH 6.2 ^m/7.7 /im depends primarily on PAH size while 
the PAH 11.3 /im/7.7/j,m is mainly sensitive to the frac- 
tion of ionized versus neutral PAHs and only modestly to 
PAH size. The effects of radiation field become noticeable 
only at high intensities (^ 10^ the average local Galactic 
far-UV flux in Habing units) and for large PAHs with 
^ 10^ carbon atoms. Mattila et al. (1999) suggested that 
PAH 6.2^m/7.7/xm variations may be due to differences 
in average temperature of the PAHs during their temper- 
ature spikes, related to PAH size or mean exciting UV 
photon energies. Alternatively, they could be attributed 
to broadening of the 7.7 /im f eatur e at low radiation field 
energy densities (Uchida et al. 200C ) ; this would reduce the 
PAH 7.7 /im flux and increase the continuum level within 
our fixed bandpasses. In addition to intrinsic variations, 
extinction can significantly alter the shape and relative 



intensities of the PAH features as shown in Sect. 5.1 



Of our sample galaxies, M82 offers the most inter- 
esting case for interpretation. The overall bilobal struc- 
ture in the maps of PAH ratios and CO gas distribution 
(Fig. ^) might result from the different composition and 
physical processes that the emitting PAHs undergo when 
exposed to the varying conditions across M82, from the 
starburst core to the more quiescent disk via the transition 
regions marked by the molecular gas ring. The maxima in 
PAH 8.6 /im/11.3 /im ratio lie at smaller radii than the 
minima in PAH 6.2 /im/7.7 /im ratio, possibly indicating 
a higher degree of PAH ionization within the most in- 
tense starburst sites at the inner edge of the molecular 
ring; this is particularly striking southwest of the nucleus 
where the peak PAH 8.6 /im/ 11.3 /im coincides very well 
with the location of the most prominent H ii region com- 



plexes (Sect. 3.2.1). The variations of PAH 6.2 /im/7.7 /im 
ratio could reflect differences in the PAH size distribu- 
tion, combined with extinction effects, where larger PAHs 
can better form and survive in denser, more shielded en- 



vironments associated with molecular gas concentrations. 
In NGC 253, the peak PAH 8.6 /im/11.3 /im ratio at the 
nucleus where the starburst is mainly occurring could also 
be due to a larger fraction of ionized PAHs. Finally, the 
enhanced PAH 11.3 /im feature in NGC 1808 compared to 
M 82 and NGC 253 could be explained by a more neutral 
mixture of PAHs and a more diffuse radiation field. 

The relative intensity of PAH and continuum emis- 
sion can also hold information on the physical environment 
within astronomical sources. In particular, the ratio of the 
peak intensity of the PAH 7.7 /im to the underlying con- 
tinuum, hereafter PAH 7.7 /im L/C ratio, has been shown 
to constitute a powerful discriminator between star for- 
mation activity or an AGN as the main source of the bulk 
of infrared luminosity (Rigopoulou et al. 1999; see also 
Genzel et al. |1998| ; Laurent et al. |200(]| ; Tran et al. ^001] ). 
Starburst-dominated objects are characterized by ratios 
^ 1 while AGN-dominated ones have ratios ^ 1 . We com- 
pared measurements of the PAH 7.7 /im L/C ratio from 
our data of M 82, NGC 253, and NGC 1808 with those ob- 



tained by Rigopoulou et al. (1999) for a large sample in- 
cluding starbursts, AGNs, and ULIRGs. Our motivation 
was to assess how much this indicator depends on the 
source luminosity at the faint end and compare the spread 
in global ratios among different galaxies with the spatial 
variations within individual galaxies. 

Figure [T^ shows the PAH 7.7/im L/C ratio versus peak 
PAH 7.7 /im luminosity for selected regions of our galaxies, 
for individual resolution elements of M 82 and NGC 253 
(within the radii defining the outer limit of the disk regions 
in table |^), and for the Rigopoulou et al. ( 1999| ) sample. 
We measured the PAH 7.7 /im peak intensity and underly- 
ing continuum according to the definitions of Rigopoulou 
et al. We note that with these definitions, the continuum 
might be underestimated in highly obscured sources as a 
result of extinction affecting notably the 11 /im region (as 
also discussed by Laurent et al. 2000). We chose the PAH 
7.7 /im peak luminosity as other characteristic property, 
taken as an approximate indicator of the total infrared lu- 
minosity. Although the PAH to infrared luminosity ratio 
can vary by up to factors of several in different environ- 
ments (Rigopoulou et al. 1999| ), our assumption has little 
consequences on the interpretation of Fig. |l2| since the 
data span a range in luminosity extending over more than 
six orders of magnitude^. 

Our data of M82, NGC 253, and NGC 1808 extend 
very well the trend defined by the global properties of 
pure starburst galaxies and ULIRGs, populating the PAH 
7.7 /im L/C ^ 1 region down to luminosities about an 



^ Additional though indirect support for our assumption 
is provided by the relationship between the 7 fim flux mea- 
sured through the ISOCAM LW2 filter (5 - 8.5 /im) and 
the far- infrare d (FIR, 43 — 123 /im) flux obtained by Roussel 
et al. (2001b) for galactic disk s, with -FV^im oc (-F'fir)^ '^. 
Furthermore, Dale et al. (2001) derived a typical value of 
« 15% for the 5 — 20 /tm to FIR luminosity ratio in normal 



galaxies while Charmandaris et al. (2002) found this ratio to 
be ~ 5% for three ULIRGs. 
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order of magnitude lower. The ratios in all three galax- 
ies lie well above the starburst-AGN separation at a ra- 
tio of unity and form a tighter distribution. The average 
and la dispersion for the individual resolution elements 
in M 82 and NGC 253 are 5.0 ± 0.9 compared to 3.7 ± 2.5 
for the pure starbursts of Rigopoulou et al. and 2.8 ± 1.0 
for their ULIRGs (excluding those with ratios < 1 or with 
only limits on the measurements). It is also interesting to 
note that the starburst trend holds for regions on spatial 
scales ranging from « 60 pc for the smallest individual re- 
gions in NGC 253 up to several kiloparsecs for the largest 
starbursts and ULIRGs (see the near-infrared images of 
Rigopoulou et al. 1999). 



6. Summary and conclusions 

We have presented ISOCAM A ~ 5.0 — 16.5 fj,m spec- 
trophotometric imaging of the starburst galaxies M82, 
NGC 253, and NGC 1808. The spectrum of all three ob- 
jects, down to the smallest scales of ~ 100 pc accessi- 
ble from our data, exhibit similar characteristics including 
prominent PAH bands, a featureless continuum most ob- 
vious at A 11 /im, and a trough in the 10 /im region. We 
securely identified the main emission features detected in 
the i? - 40 ISOCAM data of M 82 and NGC 253 based 
on their R ~ 500 — 1000 SWS spectra. The comparison 
emphasizes the caution that should be exercised when in- 
terpreting low resolution data because of potential blends 
or misidentifications between emission lines from ionized 
gas and PAH features originating in PDRs, a notable ex- 
ample being [Neii] 12.81 /im and PAH 12.7 /im. Using a 
simple model combining a template PDR spectrum and a 
power-law f^(x{X — 8.5)^'^, we constructed a representa- 
tive starburst SED and explored the effects of extinction 
at MIR wavelengths. Our simulations illustrate the impor- 
tance of the assumed extinction law (e.g. the widely used 
Draine |1989 law versus the Galactic Center law of Lutz 
1999| ) and of the intrinsic PAH spectrum (especially the 
gap between the main 6 — 9 /im and 11 — 13 /im com- 
plexes) in shaping the SED of astronomical sources. This 
complicates the interpretation of PAH ratios as well as 
extinction measurements relying on the silicate dust ab- 



sorption at 9.7 /im (see also Sturm et al. 2000). 

As observed previously in a wide range of Galactic and 
extragalactic sources, the 5 — 11 /im spectrum in our galax- 
ies is nearly invariant. The relative PAH intensities ex- 
hibit nevertheless measureable and significant variations 
of 20% — 100% which may be attributed to various, pos- 
sibly interrelated effects including the intensity of the in- 
cident radiation field and the PAH size distribution, ion- 
ization, and dehydrogenation. In our sample, M 82 proba- 
bly best illustrates variations of PAH ratios due to an in- 
creased fraction of ionized PAHs within the most intense 
starburst sites and, admittedly speculatively, perhaps also 
to differences in typical PAH sizes depending on the molec- 
ular gas concentrations. The PAH 7.7/im L/C ratio in all 
three galaxies clearly lies in the range observed for pure 



starburst systems and extends the trend reported previ- 
ously by Rigopoulou et al. ( 1999 ) to lower luminosities. 

In contrast, the A ^ 11 /im region varies most among 
our sample galaxies and the ISOCAM maps show a com- 
paratively more compact 15 /im continuum distribution 
relative to the PAH emission. Strong 15 fj,m continuum in 
M 82 and NGC 253 indicates an important contribution by 
VSGs contrary to the case of NGC 1808 where the long- 
wavelength emission is much flatter presumably because 
of negligible contribution by VSGs. We find, however, 
that in all three galaxies the 15 /im continuum and [Aril] 
6.99 /im line fluxes satisfy a linear relationship. We infer 
from this that the 15 /im continuum provides a good indi- 
cator of star formation activity in starbu rsts, co mplcment- 
ing the similar results of Roussel et al. ( ^OOlbl ) for galac- 
tic disks. In a broader perspective, our galaxies fit well in 
the trend of increasing ISOCAM 15/im/7/im ratios with 
higher levels of star formation activity found among nor- 
mal disk galaxies an d star burst-powered LIRGs / ULIRGs 
(e. g. La urent et al. 
al. ^OOll) . 



2000; Roussel et al. 20011; Dale et 



The value of the ISOCAM spectrophotometric imag- 
ing of M82, NGC 253, and NGC 1808 presented in this 
paper also lies in that it complements existing MIR data 
with maps of PAH features, fine-structure lines, and con- 
tinuum components not previously imaged. The poorer 
angular resolution of ISOCAM compared to that achieved 
with large ground-based telescopes is compensated by the 
larger field of view revealing more of the large scale emis- 
sion. In that respect, we stress the small size of the MIR 
source relative to the optical extent of all three galaxies: 
the ISOCAM maps cover the central w 1.5 kpc for M82, 
0.6 kpc for NGC 253, and 5 kpc for NGC 1808 while the 
optical diameters are about 10 kpc for M 82 and 20 kpc for 
the other two. By measuring the flux density within the 
ISOCAM LWIO filter bandpass, equivalent to the IRAS 
12 fim band, we recover all of the IRAS 12 /im emission in 
the entire field of view for M 82 and NGC 1808, and about 
70% for NGC 253. This suggests that the total MIR emis- 
sion is strongly dominated by the starburst sites in the 
nuclear regions while the more quiescent star formation 
taking place in the disk at larger radii does not contribute 
much. It will be interesting to see whether the MIR cam- 
eras on board SIRTF, to be launched in 2003, will confirm 
this result or discover faint diffuse emission, especially at 
shorter (3 — 5 /im) wavelengths. 
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Fig. 3. Mid-infrared spectra of M82, NGC 253, and NGC 1808. All spectra are normalized to unit total flux density 
between 6.0 and 6.6 /im. The three panel pairs to the left show the spectra and relative formal effective uncertainties 
for selected regions: the ISOCAM field of view, the disk region, the starburst core, and the MIR peak (see labels 
in each plot); the apertures used, normalizing factors, and typical uncertainties are given in table ^. The top and 
middle panel pairs to the right show the average spectrum together with the dispersion, the full range observed, and 
the median uncertainties derived from all individual resolution elements for M 82 and NGC 253, respectively. The 
resolution elements correspond to rebinned 2x2 pixels for M 82 and 4x4 pixels for NGC 253 (i.e. 6" x 6" in both 
cases) . The bottom right panel pair compares the spectra of the starburst core regions of each galaxy and their relative 
uncertainties. The relative uncertainties are expressed as fraction of the measured flux density. 
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Fig. 4. Comparison of MIR spectra of M 82 and NGC 253 obtained with ISOCAM and SWS. The top panels show the 
full resolution SWS spectra {R ^ 1000 to 500 from short to long wavelengths), with the identification of the detected 
emission features (from Forster Schreiber et al. |200l| for M 82, and Sturm et al. |200C| for NGC 253) . For M 82, the peak 
flux densities for the [Nell] 12.81 and [Nelll] 15.56 ^m hues are 291.8 Jy and 109.1 Jy. For NGC 253, those for 
the [Aril] 6.99 /xm, [Nell] 12.81 /.tm, and [Nelll] 15.56 ^m hues are 31.2 Jy, 140.6 Jy, and 35.2 Jy. The middle panels 
compare the SWS data convolved to the resolution of ISOCAM {R ~ 35 to 45 from short to long wavelengths) with 
the ISOCAM spectra (with uncertainties) integrated over the SWS aperture, while the bottom panels show the ratio 
of the ISOCAM and smoothed SWS spectra. In the spectral range covered by ISOCAM, the SWS aperture changes 
from 14" X 20" shortwards of A w 12.0 /j,m to 14" x 27" at longer wavelengths. The long wavelength portions of the 
spectra have been scaled to match the smaller SWS apertures by factors of 0.83 for M82 and 0.9 for NGC 253, as 
determined from the continuum level measured in the ISOCAM data. 
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Fig. 5. Selected ISOCAM maps of M82. The maps are identified at the top of each panel while the physical linear 
scale and the FWHM of the PSF are indicated at the bottom left and right. The grayscale is linear, with white and 
black tones corresponding to the minimum and maximum values displayed. For the continuum and emission feature 
images, the solid contours are 10%, 25%, 40%, 55%, 70%, 85%, and 95% of the maximum fluxes. For the 5.5 and 
15.0 fim continuum maps, the peaks are 10.3 and 187 mJy arcsec"^, respectively. For the [Aril] 6.99 /um and the PAH 
6.2, 7.7, 8.6, and 11.3 maps, they are 2.05, 12.5, 28.9, 3.84, and 5.54, respectively, in units of 10""'^^ W m~^ arcsec"^. 
The contour levels for the ratio maps are as follows: from 10 to 18 in steps of 2 for the 15.0 /xm/5.5 jim continuum 
ratio, from 0.42 to 0.48 in steps of 0.02 for the PAH 6.2 /im/7.7 /xm ratio, and from 0.45 to 0.75 in steps of 0.10 for the 
PAH 8.6 /xm/11.3 /im ratio. The dotted contours in the continuum and feature maps enclose the regions where data 
values exceed 3a, and ratio maps are displayed for the area where this is satisfied by both images involved. The cross 
at relative coordinates (0", 0") indicates the location of the galaxy nucleus. 
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Fig. 6. Selected ISOCAM maps of NGC 253. The maps are identified at the top of each panel while the physical linear 
scale and the FWHM of the PSF are indicated at the bottom left and right. The grayscale is linear, with white and 
black tones corresponding to the minimum and maximum values displayed. For the continuum and emission feature 
images, the contours are 10%, 25%, 40%, 55%, 70%, 85%, and 95% of the maximum fluxes. For the 5.5 and 15.0 /xm 
continuum maps, the peaks are 13.6 and 269 mJy arcsec"^. For the [Aril] 6.99 /zm and the PAH 6.2, 7.7, 8.6, and 
11.3 /um maps, they are 1.50, 7.66, 19.1, 2.33, and 4.08, respectively, in units of 10"^^ Wm~^ arcsec"^. The contour 
levels for the ratio maps are as follows: from 0.40 to 0.55 in steps of 0.05 for the PAH 6.2 /im/7.7/im ratio, and from 
0.30 to 0.60 in steps of 0.10 for the PAH 8.6 /xm/11.3/zm ratio. The dotted contours in the continuum and feature 
maps enclose the regions where data values exceed Scr (for the 5.5 /um continuum map, the 3cr limit lies between the 
contours of 55% and 70% of maximum flux), and ratio maps are displayed for the area whore this is satisfied by both 
images involved. The cross at relative coordinates (0",0") indicates the location of the galaxy nucleus. 
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Fig. 7. Selected ISOCAM maps of NGC 1808. The maps are identified at the top of each panel while the physical 
linear scale and the FWHM of the PSF are indicated at the bottom left and right. The grayscale is linear, with 
white and black tones corresponding to the minimum and maximum values displayed. The contours correspond to 
10%, 25%, 40%, 55%, 70%, 85%, and 95% of the maximum fluxes. For the 15.0 /xm continuum map, the peak is 
18.7 mJyarcsec-2. For the [Aril] 6.99 /zm and the PAH 6.2, 7.7, 8.6, and 11.3 urn maps, they are 0.13, 1.82, 4.23, 0.71, 
and 1.58, respectively, in units of 10~^^ W m~^ arcsec"^. The dotted contours in the contimmm and feature maps 
enclose the regions where data values exceed 3a. The cross at relative coordinates (0",0") indicates the location of 
the galaxy nucleus. 
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Fig. 9. Major and minor axis light profiles for M82, NGC 253, and NGC 1808. The curves represent the variations of 
the integrated flux per unit surface area as a function of projected distance from the nucleus for the [Aril] 6.99 fim line 
(grey lines and open circles), the 13.5 — 15.0 fj,m VSG-dominated continuum (black lines and filled circles), and the 6.0 — 9.0 
PAH-dominated region (black lines and open circles). The width of the synthetic slits along each axis is 3 pixels for M82 
and NGC 1808 and 6 pixels for NGC 253 (corresponding to 144, 477, and 108 pc, respectively). The curves are normalized 
to the maximum values which are for M82, NGC 253, and NGC 1808, respectively: 0.26, 0.36, and 0.0036 mjypc'^ for the 
"PAH" major axis profiles, 0.28, 0.38, and 0.0035 mjypc"^ for the "PAH" minor axis profiles, 0.53, 1.20, and 0.0044 mjypc"^ 
for the "VSG" major axis profiles, 0.57, 1.27, and 0.0044 mjypc"^ for the "VSG" minor axis profiles, 6.32, 7.87, and 0.035 
10"^^ Wm"^pc"^ for the [Aril] major axis profiles, and 5.77, 8.30, and 0.024 10"^^ Wm"^pc"^ for the [Aril] minor axis 
profiles. The PSF profiles for the spatially-smoothed data cubes are shown as well (simple black lines; see Sect. 0). 
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Fig. 10. Extinction effects at MIR wavelengths. The top panels show the extinction law fro m Dr aine ( 1989 ), and the 
modifications to this law at A < 10 /xm appropriate for the Galactic Center line of sight (Lutz 1999 ). The middle panels 
show the variations of the SED of a template spectrum affected by various levels of extinction, for different dust and 
source geometries and extinction laws, as labeled in each plot. The template combines the spectrum of the Galactic 
PDR NGC 7023 and a power-law spectrum (shown individually in the bottom panels by the dashed and dotted lines, 
respectively). "UFS" and "MIX" stand for uniform foreground screen and mixed model. The grey shading outlines 



the range of variations for Ay^^ =0 — 50 mag and A 



MIX _ n 
V — u 



700 mag. The various curves show the results for 



selected extinction levels, corresponding to similar attenuation factors at 10 /im between the two extinction models 
(labeled in each plot). The bottom panels show non-formal fits (grey solid line) to the spectrum of the starburst core 
of M82 (black solid line), using the template PDR -I- power-law spectrum, attenuated by the best-fit extinction for 
each geometry assuming the GC extinction law and derived from radio to optical H recombination line measurements 
(see Sect. |5.1|). All spectra are normalized to the total flux between 6.0 and 6.6 /im, except for the power-law. 
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Fig. 12. Diagnostic diagram for dominant starburst versus AGN activity. The data for selected regions in M82, 
NGC 253, and NGC 1808 as well as individual resolution elements for M82 and NGC 253 (large labeled and small 
unlabeled filled symbols, respectively; see inset) are compared to the global properties for the samples of starburst 
galaxies, AGNs, and ULIRGs of Rigopoulou et al. (1999; open circles, crosses, and stars). For M82 and NGC 253, 



resolution elements within radii of 30" and 15" (outer disk annulus; table ^ are plotted, and typical formal uncertainties 
are indicated by the error bar at the top left: < 1% for ipoak PAH7.7/im and 4% for PAH 7.7 fim L/C ratio (uncertainties 
are comparable to or smaller than the symbol sizes; average and median values differ by < 1%). 
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